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Abstract. Motivated by the recent lattice data that J /if) will survive up to 1.6T C , 
we calculate the thermal width of J ftp at finite temperature in perturbative QCD. 
The inputs of the calculation are the parton quarkonium dissociation cross sections at 
■ the NLO in QCD, which were previously obtained by Song and Lee, and a gaussian 

charmonium wave function, whose size were fitted to an estimate by Wong by solving 
the schrodinger equation for charmonium in a potential extracted from the lattice at 
finite temperature. We find that the total thermal width above 1AT C becomes larger 
C*~) ' than 100 to 200 MeV, depending on the effective thermal masses of the quark and 

gluon, which we take it to vary from 600 to 400 MeV. 
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1. Introduction 

The original claim by Matsui and SatzpQ, that J/ip suppression is a signature of 
quark gluon plasma, has witnessed a number of landmark developments, every aspect 
of which has to be taken into account consistently in confronting the recent RHIC 
data[2], and in predicting results for LHC. Among these theoretical developments are 
the phenomenologically successful statistical model for J ftp production [3j E], based 
on a coalescence assumption near T c [5], and the recent lattice calculations, showing 
strong evidence that the J/ip will survive up to 1.6T C [6], HH E]- While these two 
results seem at odd with each other, it only suggests that one still needs a more detailed 
understanding of the properties of heavy quark system in the quark gluon plasma, 
especially for temperatures between T c and 1.6T C , before a consistent picture of J/xjj 
suppression in heavy ion collision is achieved. 

In this respect, an important quantity to investigate is the effective thermal width, 
and/or the effective dissociation cross section of J/ip in the quark gluon plasma. Except 
for its existence, the present lattice results are far from making definite statements on 
the thermal width for Charmonium states above T c . Hence, in this work, we will use the 
perturbative QCD approach to calculate the thermal width. So far, such calculations 
have been limited to dissociation processes by gluons to LOpJ], [121 [13], because the 
elementary J/^-parton dissociation cross section was available only to that order [TH [15]. 
Recently, two of us have performed the dissociation cross section calculation to NLO 
in QCD[16j. Here, we will implement the NLO formula, to calculate the corresponding 
thermal width of Charmonium states in the quark gluon plasma. 
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The NLO calculation of J/^-parton dissociation calculations involves collinear 
divergence. When applying this elementary cross section to dissociation by hadrons, 
the collinear divergence is cured by mass factorization, which renormalizes the divergent 
part of the cross section into the parton distribution function of the hadron. Such 
complications disappear at finite temperature, as the thermal masses of the partons 
automatically renders the divergence finite. 

2. Quarkonium hadron interaction in QCD 

Let us begin with some introduction on the propagation of heavy quarks in the QCD 
vacuum. The propagation of a heavy quark can be approximated by a perturbative quark 
propagation with a perturbative gluon insertion, which probes the non perturbative 
gluon field configuration in the QCD vacuum. Hence, the full heavy quark propagator 
is, 

tS A (q) = iS(q) + iS{q){-ig4)iS{q) ■ ••, (1) 

where, iS(q) = i/(fa — m) and m is the heavy quark mass. The description in Eq.([T]) is 
valid even for q — > 0, because m ^> A ~ Aq C d, where in the end only gauge invariant 
combination of the gauge field A will remain after taking the vacuum expectation value. 

Table 1. Physical processes involving two heavy quarks. 

q 2 Process Expansion parameter 

Photo production of open charm 

— Q 2 < QCD sum rules for heavy quark system 
m j/V> > Dissociation cross section of bound states 



A 2 

Am 2 
A 2 

4m 2 +Q 2 
A 2 

_2QCD_ 

4rn 2 —rn 2 , , 



The propagation of a system composed of a heavy quark and an antiquark can 
also be approximated by combined perturbative heavy quark propagator with gluon 
insertions. However, since there a two heavy quarks involved, based on the operator 
product expansion, the propagation can typically be written in the following form. 

n(g) = .... + f 1 dx-— — f^/^ — r-r-7rr — (G n ) ■ (2) 

vy; Jo (4m 2 - q 2 - (2x - 1) V) n 

where, F(q 2 ,x) is a function depending on the structure of the two quark system and 
{G n ) ~ A-lfcD denotes the typical gauge invariant expectation value of gluonic operator 
of dimension 2n. The integration variable x can be thought of as the momentum fraction 
carried by one of the heavy quark. Here, one notes that such perturbative expansion 
is valid when 4m 2 — q 2 ^> Aq CD . The cases where this condition is satisfied and 
perturbative QCD treatments are possible are summarized in Table 1. In the last line of 
Tabled], 4m 2 — m 2 ,,^ ps (2m+mj/^)e , where e is the binding energy of the J/ijj. In QCD 
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Figure 1. The (a) LO diagram and (b) quark induced NLO diagram. q,pi,p2,k are 
the four momentum of the J/ip, c, c and the incoming parton respectively. 



if m — > oo, the bound state becomes Coulombic and eo = m[N c g 2 / (16tt)} 2 ^> Aqcd in the 
large N c limit. Therefore, the expansion parameter becomes small and the dissociation 
cross can be calculated using perturbative QCD. 

3. Thermal width of J/ip 

The LO invariant matrix element for the J/ip dissociation cross section by gluon was first 
obtained by Peskin[14j and rederived by one of us[T7] using the Bethe-Salpeter equation. 
The NLO invariant matrix element was derived by two of us[16]. The dissociation by 
quarks comes in as part of the NLO calculation. The thermal width of J / ip is calculated 
by folding the elementary cross section with the thermal parton distribution function. 



where deg is the degeneracy of the partons, n(k ) is the thermal quark or gluon 
distribution function and a(s) is the elementary cross section, with s being the square 
of the sum of the gluon and charmonium four momentum, and v their relative velocity. 

In both the LO and NLO calculations, a(s) is proportional to \d(p(p) / dp\ 2 , where 
<p is the charmonium wave function and \p\ is related to the incoming four momentum 
that breaks the J/ip into cc. As temperature increases, n(ko) favors higher ko or \k\, 
while the wave function favors smaller \p\, because the wave function becomes shallow 
and the momentum space wave function becomes smaller, for which we take the result 
given in ref[13j. 

For the LO calculation shown in Fig 1 (a), k = €q +p 2 /(2m c ). Hence, larger k 
means larger p. That means at higher temperature, although there are more energetic 
gluons, these are inefficient in breaking the charmonium states. Moreover, if one 
introduces thermal gluon mass, the LO width becomes even smaller. 

For the NLO calculation, Fig 1 (b) and the corresponding diagram where the 
incoming quark is replaced by a gluon, have collinear divergences, which in the vacuum 




(3) 
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is cured by mass factorization[16j. At finite temperature, the problem can be avoided 
naturally by introducing thermal masses, which we take here to be from 600 MeV[18] to 
400 MeV, the smaller value giving an upper limit to the effective width. Moreover, even 
an energetic parton can radiate a gluon with small four momentum (r in Fig 1 (b)), 
which can effectively dissociate the J/ip into open charms. Fig 2 shows the effective 
width of J/ip obtained by folding the elementary NLO quark (a) and the corresponding 
gluon (b) dissociation cross section obtained in ref[16j by a massive thermal parton 
distribution. As can be seen from the figure, even with a upper (lower) limit of 600 
(400) MeV thermal mass, the sum of thermal width due to quarks and gluons becomes 
larger than 100 (200) MeV above 1.4T C , where we take T c = 170 MeV. Hence, although 
the J/ip might start forming at 1.6 T c , the effective thermal width is very large and will 
not accumulate until the system cools down further [5]. 



(a) quark induced width 
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(b) gluon induced width 



260 280 



- m o =400 MeV 
-m°=600 MeV 



180 200 



220 240 
T (MeV) 



260 280 



Figure 2. Effective thermal width due to (a) quarks (b) gluons. Squares and circles 
are the results at temperatures where the J/ip radius is calculated in ref[13 . 
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